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1 Introduction 



During the last few years, there has been substantial experimental progress on the spec- 
troscopy of the light and heavy baryons. Among the discovered many new states, some are 
described by the quark model, but the quark content of some others is under debate. Many 
new states containing a single heavy quark are experimentally observed. Practically, all 
baryons containing a single charm quark that predicted by the quark model have already 
been established. The heavy spin-1/2, A&, E&, and Qb baryons as well as the spin-3/2, 
E£ baryon containing a single 5-quark are also experimentally discovered (for a review, 
see for example [1]). Recently, the CMS Collaboration announced the observation of the 
spin-3/2, baryon containing also a single b-quark with a mass of 5945 MeV [2]. 

The quark model also predicts hadrons with two or three heavy quarks. The experimen- 
tal progress on the heavy hadron spectroscopy has stimulated the researches on the doubly 
heavy hadron physics. Up to now, only one heavy baryon with two charm quarks, namely 
the S+ baryon has been observed in experiments [3-5] conducted by the SELEX Collab- 
oration. Researches are waiting for considerable experimental progress on observations of 
the doubly heavy baryons and their properties at LHC-b. 

One of the main characteristic parameters of the doubly heavy baryons is their mass. 
The masses of these baryons are tried to be estimated in different frameworks such as quark 
model [6, 7] and MIT bag model [8]. The masses of the spin-1/2 doubly heavy baryons are 
also calculated within QCD sum rules method [9] in [10-13], and the masses of the spin-3/2 
doubly heavy baryons are studied within the same framework in [10, 11] and [14]. However, 
one can easily see that the analytical expressions of these three works are different, and 
therefore there appears a necessity for a more accurate study of the masses and residues of 
the doubly heavy baryons within the QCD sum rules method. 

The outline of the paper is as follows. In section 2, we obtain the sum rules for the 
masses and residues of the doubly heavy baryons. In section 3, we present the numerical 
analysis of the sum rules and discuss the results. We also compare the obtained results 
with those predicted via other nonperturbative approaches in this section. 

2 Mass sum rules for the doubly heavy spin— 3/2 baryons 

The QCD sum rules for the doubly heavy spin-3/2 baryons are obtained by considering 
the two-point correlator, 



where q is the four-momentum of the doubly heavy baryon and r/ M is its interpolating 
current. 

The interpolating current for spin-3/2 baryons can be constructed from the condition 
that it should be symmetric with respect to the exchange of two heavy quarks. Using this 
condition, the interpolating current for these baryons can be written as 




(1) 



1 



e abc {(q aT C^Q b )Q' c + (q aT C^Q' b )Q c + (Q aT C^Q' b )q c ] , 



(2) 
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where a, b, c are the color indices, q is the light; and Q and Q' are the two heavy quarks. 
We present the quark contents of the doubly heavy baryons in Table 1. 



baryon 


Light quark q 


Heavy quark Q 


Heavy quark Q' 


^QQ 


u or d 


b or c 




O* 
lL QQ 


s 


6 or c 




^QQ' 


u or d 


6 


c 


O* 

12 QQ' 


s 


6 


c 



Table 1: The quark content of the spin-3/2 doubly heavy baryons 

As is well known, in the QCD sum rules approach, the correlation function is calculated 
in two different manners: 

• In terms of quarks and gluons using the operator product expansion (OPE), which 
contains perturbative and nonperturbative condensate contributions, 

• In terms of hadrons (physical part). 

Equating these two representations and performing Borel transformation with respect to 
the baryon momentum square, which suppresses the higher states and continuum contribu- 
tions, we obtain the mass sum rules. Here we would like to make the following cautionary 
note. The interpolating current f]^ of doubly heavy baryons can interact not only with 
spin-3/2 baryons, but also with spin-1/2 baryons, and this unwanted contribution should 
be eliminated. Indeed, the matrix element of the interpolating current i]^ between the 
vacuum and the heavy spin-3/2 baryons is determined as, 

(0 |^| B(q, s = 1/2)) = (a i7 „ + A 2 q,) . (3) 

Multiplying both sides with 7^ and using the fact 7^77^ = 0, we immediately get 

(0|^|S(9, a = l/2)>=A 1 ^-^)«(?), 

that is, the terms proportional to 7^ and g M contain contributions coming from spin-1/2 
baryons, and hence, such structures should be removed in further discussions. 

Keeping this cautionary remark in mind, we can now proceed to calculate the physical 
part of the correlator given by Eq. (1). Saturating this correlator by the ground state 
baryons we get, 

(0 |^| B( g ))(g(g) |^|0) 
V " g 2 -m| + '"' [) 

where dots represent the higher states and continuum contributions. The matrix element 
in Eq. (4) is determined as, 

(0\Vv\B{q)) = \u li (q) ■ (5) 
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where A is the residue and w M is the Rarita Schwinger-spinor. Using Eqs. (4) and (5) 
and performing summation over spins of the Rarita-Schwinger spinor which is given by the 
relation, 



^u^q, s)u v (q,s) 



3m% + ' 



qvip. 



3m b 



(6) 



we obtain the expression of the physical part of the correlator. As has already been noted, 
the structures containing q^, 7^, q v and j v bring contributions from spin-1/2 baryons too, 
and therefore should not be considered for determination of the characteristics of spin-3/2 
baryons. Therefore, only the structure g^ u contains solely contributions of spin-3/2 baryons 



which we shall consider in further discussions. In summary, the structures 



or 



9iw 



can 



be used for determination of the mass and residues of the heavy spin-3/2 baryons. For the 
physical part of the correlator finally we get, 



iW) = 



\ 2 {j + m B ) 
q 2 — m 2 B 



9iw + --- 



(7) 



We now return our attention to the calculation of the correlator from the QCD side. This 
calculation is carried out in deep Euclidean region using the OPE. After some calculations, 
we obtain expression of the correlator for the baryons containing two different heavy quarks 
in terms of light and heavy quarks propagators as follows: 

^^abc^a'b'c' f cc&'„, caa'„, cbc' 

3 v 

_ ocfc' qac' 



n 



Qcb' Qaa' qbc' CCi' obb' oac' G ca 7, Qbb' qo 



icb' „ G aa '„ Qbc' oca' - ebb' oac' C c ^'^ Q aa ' ^ Qbc' 
,Or> Jg lv>->Q> l/i^Q — O q ^pOq 7 m Oq/ 



IP Q 



- S$ Tr 



- ST Tr 



S Q' lv S Q 7^ 



aV 



(8) 



where S = CS T C and C is the charge conjugation operator. For the baryons containing 
the same heavy quarks, we need to replace the Q' by Q in the above equation. 

It follows from Eq. (8) that, in order to calculate the correlator from the QCD side, the 
expressions of the heavy and light quark propagators are needed. Their expressions in the 
coordinate representation are given as 



S q (x) = 
Sq(x) = 



m r , 



2n 2 x A 



An 2 x 2 



m. 



m 2 Q K x {m Q ^T- 
v'~ 



m 2 Q ft 



4vr 2 



x/ 



4n 2 x 2 



K 2 (m Q V-x 2 ) , 



(9) 



where K\ and K 2 are the modified Bessel functions of the second kind. 

It should be noted here that the propagators contain also pieces proportional to the 
gluon field strength tensor. While we perform numerical analysis with these terms, we see 
that their contributions are very small, and for this reason we do not present them in Eq. 
(9). 

The correlation function for the structure (jig^ v or g^ v in QCD side can be written in 
terms of the dispersion relation as 



n,(g 2 ) = / 



ds- 



(m Q +m Ql ) 2 



Pi 

s — q 2 



(10) 
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where i — 1(2) corresponds to the structure fig^ (g^u)- The spectral density pi in Eq. (10) 
is given by the imaginary part of the correlator, 



Pi(s) = -Imrij(s) 



7T 



In calculating the correlator from QCD side, we keep the term linear in strange quark 
mass. After tedious calculations, for the spectral densities we get 

"j rCX-max r ft max s 

Pi{s) = / da dfilyp 3a(3(a + (3)p - 4(1 - a - 0)m Ql mQ 

- Am q (am Q > + j3m Q ) j - -^-^ / da (1 - a)(3am q - 4m Q ) - 4am Q/ (11) 

P2(s) = / da d(3^p^(a + [3)p(am Q , + (3m Q ) - m q (a(3p + 3m Q mQ/)J j 

*/ CKrn.ir>. ** Pm.i.m. 



+ 



(gg) 
12tt 2 



x-mtn 

/ da|(l - a) 2a(mo + //i - s) -mgm g + (am, + 3m Q )m Q/ j , (12) 



where, 



/i = 



- s , 



a /3 
/ii = /i(/3 ->■ 1 - a) , 



am; 



sa — m 



2 ' 





1 — a , 
1 



av 



a r 



2s 
1 ' 

2s 



, 2 2 

s + m -m c 



, - y / '(s + m 2 Q - m 2 Q ,) 2 - Am 2 Q s , 



s + — m%, + \/(s + m?Q — m 2 n ,) 2 — 4m%s 



We can now compare our results on the spectral densities with the ones presented for 
instance in [10] and [11]. As far as the spectral density pi(s) is concerned, we have the 
factor 3(a + j3), which is absent in [11]. The quark condensate in our case contains the 
term proportional to a ^~^ m i ; while the corresponding term in [11] is ||ai(l — a ) m q- The 
differences for the spectral density p2(s) can be summarized as follows. The perturbative 
term proportional to mq{mQi) (note that rnqimQi) in our work corresponds to the m,Q'(mQ) 
in [11]) contains the factor (a + f3) which is absent in [11]. In the perturbative part of p2(s) 
we have also 



1 [da [dp. 2 2 2 



which is again absent in [11]. When we compare the quark condensate terms we have 
4cc(l — a)s which is different from their term reading 3a(l — a)s. 

structure, our results on perturbative part and quark condensate terms 



For the 



without the strange quark mass agree with the ones given in [10], but the terms proportional 



4 



to m s (which is calculated in [15]) and the results for the d — 5 operators are different 
compared to those given in [10] and [15]. 

Equating the coefficient of structure ig^Qnu) in Eq. (7) to that of Eq. (10) for Ili(n 2 ), 
and performing Borel transformation with respect to Q 2 = —q 2 , we get the sum rules for 
the mass, 



J(m Q +m' Q ) 2 

in which the quark-hadron duality is used, and the contributions of the higher states and 
continuum are modeled as the perturbative ones starting from some threshold s . 



In this section, we present our numerical results on the mass and residues of the spin- 
3/2 doubly heavy baryons. For the quark masses, we use their MS values: m c (m c ) = 
(1.28 ± 0.03) GeV, rh b (m b ) = (4.16 ± 0.03) GeV (see for example [16]), and m s {2 GeV) = 
(102 ± 8) GeV [17]. The values of the quark condensates are taken as (uu)(l GeV) = 
(dd)(l GeV) = -(246+^ MeVf [18], (ss) = 0.8(uu) and m 2 = (0.8 ± 0.2) GeV 2 . 

According to the sum rules analysis, the working regions of the continuum threshold Sq 
and the Borel mass M 2 should be found by imposing the requirement that the mass and 
residue exhibit good stability with respect to the variations in these parameters. Therefore, 
we vary the continuum threshold So and the Borel mass M 2 , in order to find the "working 
region" of M 2 , where the perturbative contribution is larger compared to the nonpertur- 
bative part. Using the quark-hadron duality, the contributions of the higher states and 
continuum are taken as the perturbative ones starting from sq. 

The continuum threshold depends on the energy in the vicinity of the first excited 
state. In this respect, we choose the value of the continuum threshold within the interval 



s = (100 -T- 125) GeV 2 for bb, s = (50 -r- 65) GeV 2 for be, and s = (14 + 22) GeV 2 for cc 



In the analysis of QCD sum rules, two conditions are satisfied for M 2 . a) The pole 
dominance with respect to the higher states and continuum; b) The convergence of the op- 
erator product expansion, i.e., dominance of the perturbative part over the nonperturbative 
contributions. 

The upper bound on M 2 can be obtained from the condition (a). For this purpose, we 
introduce the ratio R, which describes relative contributions of the continuum and pole, 



J \m Q +m' Q ) 2 

Demanding that R < 1/2, which guarantees that the pole contribution exceeds the contin- 
uum and higher state contributions, we find the maximum values of M 2 for cc, be and bb 




(13) 



(14) 



3 Numerical analysis 



baryons. 
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baryons, as are listed below, 



( 4.5 GeV 2 (at ^ = 4.4 GeV), for 2* c and Q* c 
M 2 ma , = { 8.0 GeV 2 (at ^ = 8.0 GeV), for ~* c and {l* bc (15) 
[ 12.0 GeV 2 (at ^ = 10.9 GeV), for ~* bb and ft£ 6 . 

The lower limit of M 2 can be obtained when the criteria (b) is satisfied. Our numerical 
analysis leads to the following minimum values of M 2 , 

( 3.0 GeV 2 (at ^ = 4.4 GeV), for E* cc and Q* cc 
M^n = { 6.0 GeV 2 (at ^ = 8.0 GeV), for S 6 * c and fi 6 * c (16) 
{ 8.0 GeV 2 (at ^ = 10.9 GeV), for E* bb and Q£ b . 

In these ranges of the M 2 , the relative contributions of the pole and continuum for each 
baryon are presented in Table 2. 





Pole 


Continuum 


1 cc 


(63 


- 71)% 


(29- 


37)% 


"cc 


(75 


- 81)% 


(19- 


25)% 


"* 

"6c 


(68 


- 75)% 


(25- 


32)% 


0* 

"&c 


(77 


- 83)% 


(17- 


23)% 


"* 
"66 


(52 


- 58)% 


(42- 


48)% 


O* 
"66 


(66 


- 70)% 


(30- 


34)% 



Table 2: The relative contributions of the pole and continuum to the sum rule in respect 
to the variation of M 2 in its "working region" . 

Using the the working regions of M 2 and So, we obtain the results for the masses of 
spin-3/2 doubly heavy baryons, which are all presented in Table 3. For completeness, we 
present the results of the other works in the same Table as well. The residues of these 
baryons are also presented in Table 4. We see from these Tables that our results on the 
masses are overall very close to the values given in [6], [10], [11], [14] and [15]. For the 
masses of and Ql b , the predictions of all approaches are very close to each other. On 
the other hand, for the masses of and Ql c , the predictions of [11] are slightly larger 
in magnitude while the results of [6] are slightly smaller compared to our results. As far 
as the residues of the E* c and fl* c baryons are concerned, our predictions are higher in 
magnitude compared to the ones presented in [10] and [14]. Our prediction for the residue 
of the baryon is also higher when compared to that of the [10], while our predictions 
on the residues of the S£ fe and Ql b baryons almost match with those of the [14] . 

In conclusion, the masses and residues of the spin-3/2 doubly heavy baryons have 
been calculated within the QCD sum rules method. A comparison of our results with 
the predictions existing in the literature has also been presented. We expect that these 
predictions can all be checked at LHC in near future. 
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Our Work 
Structure $g^ u Structure g^ v 


[10] and [15] 


[11] 


[6] 


[14] 


1—1 cc 


3.o9±0.1o 


3.7z±U.18 


3.58±0.05 


3.90±0.10 


3.727 


o a i i o 1 o 

3.ol±0.18 




3.78±0.16 


3.78±0.16 


3.67±0.05 


3.81±0.06 


3.872 


3.76±0.17 


^bb 


10.4±1.0 


10.3±0.2 


i n qq+i f)Q 


in Qtjxn no 


1 D 237 


iri 99+0 1 5 


^cb 


7.25±0.2 


7.2±0.2 




8.00±0.26 


6.98 




o* 


7.3±0.2 


7.35±0.2 




7.54±0.08 


7.13 




0* 


10.5±0.2 


10.4±0.2 


10.38±1.10 


10.28±0.05 


10.389 


10.38±0.14 



Table 3: The mass spectra of the spin-3/2 doubly heavy baryons (masses are in units of 
GeV) 





Our Work 
Structure ^<7 Mi , Structure g^ u 


[10] 


[14] 




0.12±0.01 


0.12±0.01 




-7.8) x lO- 3 


0.070±0.017 


0* 


0.14±0.02 


0.13±0.01 






0.085±0.019 


"* 
"66 


0.22±0.03 


0.21±0.01 


V(5-0- 


23.0) x 10- 3 


0.161±0.041 


"* 


0.15±0.01 


0.15±0.01 








o* 


0.18±0.02 


0.17±0.01 








O* 
s '66 


0.25±0.03 


0.25±0.02 






0.199±0.048 



Table 4: The residues of the spin-3/2 doubly heavy baryons in units of GeV 3 . 
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